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ABSTRACT 


The  operation  of  a  bolometer  thermal  sensor  is  analogous  to  that  of  a  charging 
capacitor  in  a  generic  RC  circuit.  As  such,  circuits  containing  bolometers  can  be 
analyzed  with  standard  circuit  simulation  programs  such  as  PSPICE.  This  thesis  deals 
with  the  development  of  a  bolometer  model  by  using  PSPICE  with  the  aid  of  Analogue 
Behavior  Modeling  (ABM)  capability,  which  allows  the  user  to  program  circuit 
components  with  basic  maAematical  functions.  The  predictions  of  the  model  were 
found  to  be  in  good  agreement  with  the  reported  experimental  data,  which  demonstrates 
the  accuracy  of  the  model.  The  model  was  used  to  design  a  self-heating  compensated 
thermal  sensor  with  enhanced  signal  integration  capability  to  improve  the  signal-to-noise 
ratio.  We  believe  the  model  can  be  used  to  analyze  any  circuit  containing  bolometers  to 
optimize  the  performance. 
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I.  INTRODUCTION 


Infrared  (IR)  technology  remains  vital  to  military,  law  enforcement  and  fire¬ 
fighting  applications.  Numerous  means  of  detecting  IR  radiation  exist,  each  with  its  own 
drawbacks  and  advantages.  Among  these,  for  example,  the  HgCdTe  and  InSb-based 
photodetectors  typically  have  superior  performance  and  moderate  cost,  but  they,  as  well 
as  many  others,  require  cooling  to  reduce  dark  current  (Rogalski,  1994).  This  in  turn 
requires  a  larger,  heavier  overall  package  and  higher  power  consumption. 

The  need  for  small,  lightweight,  uncooled,  sensitive  IR  devices  has  grown  in 
particular  in  the  military  as  a  complement  to  the  trend  of  smaller  unit  tactics  and  in 
response  to  the  need  for  more  intelligence  and  identification  capability  at  the  organic 
level.  From  imaging  devices  motmted  on  rifles  to  thermal  detectors  disguised  as  small 
plants  providing  perimeter  surveillance,  uncooled  IR  technology  holds  the  promise  of 
providing  an  edge  over  the  enemy. 

The  thermistor  bolometer,  which  we  will  refer  to  simply  as  the  bolometer,  is  a 
temperature-variable  resistor  which  can  be  used  in  imcooled  systems.  By  measuring  the 
change  in  resistance,  a  signal  can  be  found  which  is  dependent  only  on  the  change  in 
temperature  which  the  bolometer  experiences  due  to  incident  IR  flux.  A  number  of 
difficulties  have  been  overcome,  such  as  vacuum-encapsulating  the  bolometer  to 
eliminate  unwanted  conduction  of  heat  from  the  sensor  by  surrounding  air.  Nevertheless, 
a  large  obstacle  to  high  sensitivity  remains:  the  signal  itself  is  affected  by  the  very  act  of 
measuring  it. 

To  measure  the  change  in  resistance  which  has  been  caxised  by  heating  the 

bolometer,  a  bias  voltage  must  be  applied  across  it.  The  resulting  current  may  then  be 
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measured,  and  the  change  in  resistance  may  be  calculated.  The  difficulty  in  measuring 
the  optically-induced  change  in  resistance  comes  about  because  of  the  additional  ohmic 
heating  that  the  current  generates.  If  we  calculate  the  incident  power  due  to  a  typical  IR 

signal  (~  nW)  versus  this  self-heating  (~mW),  also  called  Joule  heating,  we  find  that  at 

detector  saturation  the  Joule  heating  is  typically  three  orders  of  magnitude  higher  than  the 
IR  signal. 

The  issue  then  becomes  how  to  discriminate  the  signal  from  the  self-heating 
response.  In  most  applications,  the  self-heating  response  is  first  measured  in  the  absence 
of  an  IR  signal.  This  dark  current  value  is  then  subtracted  fi-om  subsequent  measured 
signals.  While  this  technique  is  relatively  inexpensive  and  easily  implemented,  it  is 
prone  to  error*  and  limits  the  dynamic  range  of  the  amplifying  circuitry. 

A  more  promising  technique  is  the  use  of  a  reference  bolometer,  which,  if 
incorporated  into  a  Wheatstone  bridge  arrangement  (see  figure  II.7),  can  theoretically 
cancel  out  the  self-heating  response  if  it  is  optically  isolated  fi-om  incident  IR  radiation. 
One  method  was  introduced  in  a  European  patent  (Kimura,  et  al,  1997),  which  provides  a 
physical  shield  for  the  reference  bolometer.  This  technique  may  be  effective,  but  it  is 
very  costly  when  implemented  into  micro-machined  bolometer  circuits. 

The  necessity  for  low-cost,  easily  fabricated  silicon  chip  technology  prompted 
Karunasiri,  et  al,  to  investigate  the  possibility  of  rendering  the  reference  bolometer 
translucent  to  IR  by  means  of  adjusting  its  thermal  characteristics.  Promising  physical 
experiments  validated  this  concept  (Karunasiri,  1999),  but  it  has  not  been  fully  explored. 

*  an  average  value  is  subtracted  from  the  actual  integrated  signal  value,  for  example. 
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Analytical  models  of  the  process  can  be  derived  using  first-order  approximations,  but 
what  was  needed  next  was  a  suitable  computer  model  that  could  account  for  higher-order 
effects  in  order  to  further  refine  the  concept  of  manipulating  the  reference  bolometer’s 
thermal  characteristics  to  improve  self-heating  cancellation. 

The  purpose  of  this  thesis  is  to  develop  a  tool  to  conduct  parametric  studies  of 
bolometer  operation.  What  follows  is  the  development  of  such  a  model  using  PSPICE, 
an  industry-standard  electrical  engineering  circuit  analysis  tool.  PSPICE  is  used  in  a 
novel  way  to  simulate  the  thermal  characteristics  of  the  bolometer.  The  model  is  verified 
against  physical  experimental  data,  and  then  an  example  is  given  of  how  the  model  can 
be  used  to  explore  and  refine  the  concept  of  self-heating  cancellation  through  the  use  of  a 
compensating  bolometer. 
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II.  THEORY  OF  BOLOMETER  OPERATION 


A.  GENERAL  MODEL  OF  THERMAL  DETECTORS 
1.  The  Heat  Transfer  Equation 

The  operation  of  a  generalized  thermal  detector  can  be  understood  by  using  a 
coimnonly  accepted  model  wherein  the  sensor  element  is  connected  to  a  heat  sink. 
Energy  absorbed  by  the  sensor  is  transported  via  a  bridge  with  thermal  conductance  G  to 
the  heat  sink.  From  the  heat  sink,  the  energy  is  either  removed  by  a  thermo-electrical 
(TE)  cooling  process  or  it  dissipates  as  the  heat  sink  re-establishes  thermal  equilibrium 
with  its  environment. 


Sensor  Heat  Sink 


A  natural  trade-off  exists  between  the  signal  response  level  of  tiie  sensor  and  how 
quickly  the  sensor  can  acquire  a  new  signal;  both  are  dependent  on  characteristics  of  the 
sensor  material  as  well  as  the  physical  connection  to  the  heat  sink.  The  process  described 
here  is  governed  by  the  heat  transfer  equation,  which  is  given  by  (Boreman,  Dereniak, 
1996) 

P^‘GAT  +  Hi^.  (II.1) 
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The  heat  transfer  equation,  which  is  a  statement  of  energy  conservation,  states  that  the 
incident  power  (Pmcuient)  equals  the  thermal  conductance  (G)  times  the  change  in 
temperature,  plus  the  heat  capacity  (H)  times  the  rate  of  change  of  temperature.  Thermal 
conductance,  not  surprisingly,  is  the  ability  of  a  device  to  conduct  thermal  energy  and  has 
units  of  watts  per  kelvin  (W/K).  Heat  capacity  (H)  is  the  ability  of  a  device  to  accept 
thermal  energy  in  the  form  of  heat  and  has  units  of  joules  per  kelvin  (J/K).  Typical 
values  of  G  and  H  for  a  microbolometer  are  10'^  W/K  and  10’^  J/K,  respectively. 


2.  Steady-state  Operation 


Assuming  a  scenario  wherein  the  sensor  is  subject  to  a  constant  incident  radiative 
power  flux,  the  sensor  will  reach  a  certain  maximum  temperature  as  its  heat  capacity  is 
filled  up.  At  this  point  the  incident  power  simply  passes  through  the  sensor  to  the  heat 
sink  where  it  is  dissipated.  This  change  in  temperature  can  be  extracted  from  the  heat 
transfer  equation  by  setting  the  rate  of  temperature  change  equal  to  zero  and  solving  for 
AT. 


When  dAT/dt  =  0, 


(11.2) 


3.  Transient  Operation 

dlsJT 

Given  the  governing  heat  transfer  differential  eqiiation  H - =  P-G-AT  ,  we 

dt 

wish  to  solve  for  the  function  AT  in  order  to  understand  the  bolometer’s  change  of 
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temperature  as  a  flmction  of  time  during  the  transient  stages.  Comparing  the  heat  transfer 
equation  to  the  equation  governing  the  charging  of  a  capacitor  in  a  simple  RC  circuit  with 
a  constant  cxurent  source  I,  note  the  similarity: 

H  —  +  G-AT  =  P  versus  C-^+  — =  /.  (11-3) 

dt  dt  R 

The  RC  circuit  which  equation  (II.3)  describes  is  shown  below. 


Figure  11.2  RC  parallel  circuit 

The  solution  to  the  RC  circuit,  assuming  the  initial  value  F=0  at  t=0,  is  F  =  IR(\-e 
and  by  comparison  we  should  expect  a  similar  solution  to  the  heat  transfer  equation  of  the 

form  AT  =  — (1  - e  if  AT=0  at  t=0: 

G 

Rearranging  the  heat  transfer  equation, 

dM  ^  -dt 
GAT-P ~  H  ' 


Lety  =  GAT -P  and  dy  =  GdAT .  Then, 


y  H 


Moving  G  to  the  right  side  and  integrating  both  sides. 
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Assigning  =  r ,  the  time  constant, 

ln(;;)  =  -— + c ,  where  c  is  an  arbitrary  constant. 

T 

To  solve  for  the  constant  of  integration,  we  apply  the  boimdary  condition  that  at 
r  =  0  the  change  in  temperature  (AT)  also  equals  zero,  and  hence  y  =  -P  (from 


y  =  GAT-P): 


In(-P)  =  c,  and  therefore  ln(y )  =  - — + ln(-P) . 

T 


Collecting  terms. 


In 


\-p^ 


and  solving  fory  yields 


y  =  -Pe  . 

Substituting  the  original  value  of  y  =  GAT-P  we  arrive  at 

GAT-P  =  -Pe'^\ 


and  finally 


AT 


where  t  =  ^/ 


G‘ 


(11.4) 


Typical  values  of  x  are  on  the  order  of  10'^  seconds'’  and  are  obviously  material  and 
design  specific. 


If  the  incident  radiation  is  removed  after  the  sensor  reaches  saturation,  P  =  0,  and 
we  may  observe  how  the  device  cools: 

0  =  GAT  +  H(dAT/dt). 

Rearranging, 

dM  ^  -dt 
GAT~  H  ’ 


and 


substituting  y=Gidr  this  time. 


where  P  is  the  power  used  to  achieve  the  steady  state  from  which  the  temperature 
then  falls.  The  following  figure  illustrates  the  relationship  between  the  change  in 
temperature  and  the  application  of  constant  power  which  is  then  turned  off 
instantaneously  (step  response). 


Elapsed  Time 


Figure  II.3  General  thermal  detector  step  response. 
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4. 


Practical  Considerations  and  Trade-Offs 


As  equation  II.  1  shows,  a  bolometer  measures  incident  power  by  changing  its 
temperature.  The  physical  parameter  that  is  actually  measured  is  the  temperature  change 
AT.  A  value  for  P  may  then  be  extracted  by  applying  equation  (II.2),  assuming  G  and  H 
are  known  quantities.  The  other  assumption  made  here  is  that  we  are  measiuing  AT  once 
the  sensor  has  reached  its  steady  state.  Otherwise,  as  figure  II.2  illustrates,  we  measure  a 
time-dependent  transient  response  akin  to  equation  II.4,  which  is  a  more  complicated 
function  of  P. 


To  describe  how  the  sensor  responds  to  incident  IR,  we  define  the  temperature 
responsivity  of  the  detector  by 

AT 


Rr  =■ 


(II.6) 


Noting  that  AT  =  P/G  for  steady-state  conditions  [equation  (II.2)],  we  may  re-write 
equation  (II.6)  as 


Rj.=  — 

G 

and  we  see  that  responsivity  is  inversely  proportional  to  G. 


(II.7) 


Another  important  parameter  for  a  detector  is  speed,  or  temporal  frequency 
response.  In  the  derivation  of  equation  (II.4),  we  determined  that  the  time  constant 

7  =  ^/^,  where  and  4/^  is  the  effective  frequency  response  bandwidth. 

Considering  this  and  equation  (II.  7)  we  see  that  speed  and  responsivity  are  inversely 
proportional.  One  gives  up  frequency  response  for  greater  responsivity  in  thermal 
detectors.  Indeed,  it  is  evident  that  the  responsivity-bandwidth  product  RjAfis  a  constant 
for  a  bolometer  of  a  given  heat  capacity. 


10 


B.  INTRODUCTION  TO  THE  BOLOMETER 

1.  Description 

A  bolometer  is  a  device  whose  electrical  resistance  changes  as  a  function  of  its 
temperature.  Bolometers  are  also  known  as  thermistor  bolometers,  or  thermally  sensitive 
resistors.  The  bolometer’s  ability  to  change  resistance  is  described  by  its  temperature 
coefficient  of  resistance  (TCR),  which  will  be  hereafter  denoted  by  a.  Typical  a  values 
for  metals  material  run  in  the  region  of  approximately  0.3%  change  per  degree  Celsius, 
but  semiconductors  may  have  values  as  high  as  2%  change  per  degree  Celsius.  The  value 
of  a  is  used  to  determine  resistance  of  a  bolometer  at  a  given  temperature  as: 

i?(r)  =  i?o(l  +  a-[7'-ro]),  (II.8) 

where  Rq  is  the  bolometer’s  resistance  at  ambient  conditions.  To  is  the  ambient 
temperature  and  T  is  the  instantaneous  bolometer  temperature.  (T-To  will  be  designated 
by  AT  hereafter.) 

A  microbolometer  is  simply  a  miniature  bolometer  micro-machined  or  photo- 
etched  onto  a  substrate,  typically  silicon.  Figure  II.4  shows  a  pair  of  microbolometer 
images. 
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Figure  II.4  Photographs  of  a  low  thermal  conductance  microbolometer  (left)  and  a 
high  thermal  conductance  microbolometer  (right).  (Ramakrishna,  et  al,  1999) 

The  bolometers  shown  above  are  identical  in  every  way  with  the  notable  exception  being 
the  addition  of  two  wide  connections  to  the  substrate  heat  sink  on  the  right  bolometer. 
These  additional  bridges  provide  more  material  though  which  heat  may  flow  to  the  heat 
sink.  Hence,  the  right  bolometer  has  a  higher  thermal  conductance. 

2.  Bolometer  Design 

From  both  engineering  and  economical  points  of  view,  it  is  of  great  importance 
that  the  entire  sensor  be  designed  to  fit  onto  a  single  silicon  chip  as  an  integrated  circuit. 
Attempting  to  attach  leads  to  individual  circuit  elements  will  otherwise  dramatically 
reduce  the  manufacturing  yield  and  increase  the  cost  per  detector.  A  cursory  look  at  the 
micromachining  process  of  a  bolometer  may  provide  a  more  thorough  understanding  of 
its  physical  characteristics.  Figure  II.5  (below)  depicts  a  typical  design. 


50  um 


50  |j.m 


Anisotropic 

Etching 


Figure  II.5  Microbolometer  pixel  design. 


The  surface  features  are  etched  using  standard  lithographic  processes.  In  order  to 
thermally  isolate  the  bolometer  (with  the  exception  of  its  bridges  to  the  heat  sink),  the 
region  below  the  bolometer  is  etched,  for  example,  using  tetramethyl  ammonium 
hydroxide  (TMAH)  solution.  This  solution  interacts  with  the  p-type  (100)  silicon 
substrate  in  a  useful  fashion:  it  only  etches  the  silicon  in  the  1 1 1  crystal  plane  direction. 
These  111  planes  can  be  clearly  seen  in  the  following  scanning  electron  microscope 
photograph  (figure  II.6).  During  measurement  the  device  is  placed  under  vacuum  to 
eliminate  any  convective  heat  transfer.  The  importance  of  vacuum  encapsulating  the 
bolometer  is  demonstrated  in  figure  II.7,  which  is  a  series  of  images  taken  by  a  focal 
plane  array  of  bolometer  pixels  as  it  is  placed  under  increasing  vacuum.  The  image 
quality  is  improved  because  the  bolometers  suffer  from  increasingly  less  convection  heat 
loss  under  increasing  vacuum  conditions. 
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Figure  II.6  Scanning  electron  microscope  (SEM)  image  of  a  microbolometer  showing 
the  etches  along  the  1 1 1  plane.  (Gu,  et  al,  1998) 


Figure  IL7  Images  taken  by  a  microbolometer  FPA  under  increasing  vacuum  of 
values  1)  ambient  pressure  2)  100  torr  3)  5  torr  4)  1x10^  torr.  (He,  et  al,  2000) 


C.  READ-OUT  CIRCUIT  USING  A  WHEATSTONE  BRIDGE 

1.  Purpose 

The  Wheatstone  bridge  (Dereniak,  1996)  is  commonly  used  for  isolating  circuit 
components  and  making  measurements  on  them.  As  a  first  step  in  employing  a 
bolometer  to  measure  temperature  or  infrared  radiation  (IR),  the  Wheatstone  bridge 
provides  a  convenient  circuit  architecture.  As  will  be  discussed  in  section  C.7,  using  a 
second  bolometer  for  compensation  improves  the  output  signal,  but  for  now  it  is 
instructive  to  look  at  a  Wheatstone  bridge  with  a  single  bolometer  element  in  it,  as 
depicted  below. 


Figure  II.  8  Wheatstone  bridge  with  one  bolometer.  In  the  text,  Vsias  will  be  denoted 

by  Vb  and  AVout  will  be  denoted  by  Fout 

2.  Signal  Voltage  as  a  Function  of  Temperature  Change 
As  previously  stated,  the  bolometer’s  resistance  changes  as  a  function  of 
temperature  in  the  following  manner:  R(T)  =  Ro(l+a-AT).  Provided  that  no  current  flows 
through  the  nodes  connecting  Fj  and  F^,  we  may  consider  both  sides  of  the  Wheatstone 
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bridge  as  voltage  dividers  with  a  bias  voltage  (F*)  applied,  in  which  case  the  voltages 
with  respect  to  ground  may  be  calculated  as  follows. 


V  =-^F 


^2  = 


R(T) 


R,+R(T) 


(II.9) 


In  an  actual  circuit  we  may  measure  the  voltage  across  the  bridge  (Vom),  which  can  be 
calculated  as 


V  =V -V =V 

^  out  ^  b 


i?Q  +  /?g  (l  +  ot^T  ^  2 


V  =F 

'out  'b 


l+aAT  1 

2+aAr  2 


V  =v 

'out  'b 


l  +  oAT 


2+Q^r 


2-aM 

l-aST 


(11.10) 


V  =v 

'out  'b 


2+a^T-a^Ar 

4-a^AT^ 


1 

2 


where  the  quadratic  terms  may  be  considered  negligible  if  \aAT\ «:  1 .  Then, 


Vout=V, 


2  +  aAT 


and  finally  we  are  left  with  a  relatively  simple  expression  for  the  signal  voltage  (Voui)  in 
terms  of  the  change  in  temperature: 


(11.11) 
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From  equation  (II.5)  we  understand  that  AT  is  an  exponential  in  time  with  a  linear 
dependence  on  incident  power,  which  means  that  in  the  Wheatstone  bridge,  Vout  is  also  an 
exponential  with  a  linear  dependence  on  power.  Two  forms  of  power  influence  the 
temperature  of  the  bolometer.  Incident  power  is  due  to  the  incident  IR  radiation,  which  is 
what  we  wish  to  measure.  The  other  form  is  dissipated  power  due  to  the  bias  voltage, 
which  is  applied  in  order  to  pass  current  through  the  bolometer  in  order  to  measure  its 
changed  resistance. 

3.  Power  dissipated  across  bolometer 

We  now  wish  to  consider  more  carefully  how  the  heating  of  the  bolometer  occurs 
due  to  an  external  bias  voltage  (Vbias).  This  is  known  as  Joule  heating  or  self  heating.  To 
calculate  the  quiescent  power  dissipated  across  the  bolometer,  we  square  the  voltage 
across  the  bolometer  and  divide  by  its  resistance: 

Pj  ,  — ,  where  F,  =  — and  R{T)  =  Rf,(l+aAT) . 

Joule  2 

P  ^(3^)  p2  _ _ m+ccAT) _ 

{Ro+R(T)f  '  R^+2R,(R^(l  +  aAT))  +  {R,(l  +  ocAT)f  ' 

P  R^jl+ocAT)  2 

4R^+4R^aAT  +  R^a^AT^  ” 


Dividing  through  by  Ro  and  treating  the  quadratic  term  as  negligible  for  oAT  4 , 


F/  l+aAT 
4Ro*l  +  aAr’ 


(11.12) 
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and  we  see  that  the  power  dissipation  due  to  Joule  heating  is  independent  of  AT  as  long 
as  |«Ar|  <g:  4 ; 


P  =-^ 

Joule 


4R, 


(11.13) 


Once  steady  state  has  been  reached,  all  of  the  incident  power  will  move  through 
the  bolometer  to  the  heat  sink,  as  can  be  described  by  the  following  set  of  equations: 


j  A  7^  j  A  7^ 

H - =  Pi„cide„,  -  gat  and  in  a  steady  state, - =  0 ,  so 

at  dt 

P,.^‘GAT.  01.14) 

As  previously  stated,  equation  (11.13)  expresses  power  dissipated  by  the 
bolometer  due  only  to  the  bias  voltage.  In  an  application  setting  in  which  we  wish  to 
measure  the  incident  IR  radiation,  we  must  recognize  that  the  voltage  response  of  the 
Wheatstone  bridge  will  be  due  to  the  sum  of  the  incident  IR  power  and  the  dissipated 
power  due  to  Joule  heating. 

4.  Wheatstone  Bridge  Response  to  Bias  Voltage 

Before  determining  the  effects  of  external  radiation,  it  is  instructive  to  look  at  the 
output  voltage  of  the  circuit  as  a  function  of  time,  due  to  an  applied  bias  voltage  alone. 
Figure  II.8.  Given  equations  (II. II)  and  (II.5),  we  may  now  write 

where  t  is  the  time  elapsed  since  the  application  of  Vb.  In  the  absence  of  incident 
radiation,  total  power  dissipated  by  the  bolometer  is  given  by  the  equation 
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(11.16) 


P  = 

Joule 


and  therefore  the  output  voltage  of  the  Wheatstone  bridge  is 


16-i2o-Gl  / 


(11.17) 


The  response,  not  surprisingly,  takes  on  the  same  graphical  form  as  that  of  the 
exponential  rise  in  temperature  depicted  in  figure  II.3.  Let  us  look  more  closely  by 
inserting  some  typical  values  into  equation  (11.17).  Figure  II.8  represents  the  output 
voltage  of  just  such  a  circuit  (parameter  values  are  given  in  the  figure  caption). 


Time  (ms) 


Figure  II.9  Output  voltage  for  a  one-bolometer  Wheatstone  bridge  with 

G  =  3.2x10'^  W/K,  H  =  4x\0'‘J/K ,  i?  =  16A:Q,anda=0.27%. 

For  t »  r  (steady  state)  the  output  voltage  asymptotically  approached  the  value  of  the 
prefix  term  in  equation  (11.17),  which  is  called  the  saturation  voltage 


saturation 


V^-cc 

\6-R^-G 


(11.18) 
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and  in  the  example  of  figure  11.9  is  calculated  to  be  approximately  10.4  mV.  Equation 
(11.18)  may  be  used  to  determine  the  bolometer  detector  system’s  G  firom  a  measurement 

of  ^saturation- 


It  is  also  useful  to  determine  an  empirical  value  for  the  system’s  //parameter.  For 
t «  T,  we  expand  the  exponential  term  to  the  first  order  and  simplify: 


\6R^G 


V!'cc  t 

16/?o-Gr’ 


and  using  equation  II.4, 


V  = 

'  out 


V,-cc 


■t. 


(11.19) 


From  this  equation  we  observe  a  linear  response  for  small  t  which  is  independent  of  G.  If 
we  were  to  collect  a  number  of  data  points  for  this  time  region,  we  would  obtain  a  plot  of 
the  form  found  in  figure  II.  10. 


Figure  II.  10  Depiction  of  the  output  voltage  response  for  Wheatstone  bridge  over  a 

period  of  time  / r . 
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The  fitted  slope  of  the  corresponds  to  the  coefficient  of  t  in  equation  (11.19),  which  can 
then  be  used  to  solve  for  H.  Hence,  one  can  determine  both  H  and  G  with  one  set  of 
measurements. 

5.  Radiative  Thermal  Conductivity 

One  important  consideration  should  be  noted,  however,  concerning 
thermal  conductance.  Heat  transfer  occurs  in  three  forms:  convection,  conduction  and 
radiation.  With  the  bolometer  imder  vacumn,  convection  is  negligible.  Nevertheless, 
while  heat  put  into  the  bolometer  will  be  conducted  through  the  bridges  attaching  it  to  the 
heat  sink,  it  may  also  radiate.  When  an  object  heats  up,  it  radiates  in  accordance  with  the 
Stefan-Boltzman  Law  (Dereniak,  1996),  which  gives  the  radiant  power  as 

=  A€C7(r -T^) ,  (11.20) 

where  To  is  the  ambient  temperature,  A  is  the  radiant  surface  area,  emissivity  s  is  a 

fimction  of  the  bolometer  material,  and  the  Stefan-Boltzman  constant 
cr  =  5.67xlO‘*w<3«/(/M^-A’'‘).  A  true  blackbody  has  an  emissivity  of  exactly  one. 
Considering  that  most  metals  are  more  reflective  than  emissive,  one  would  assume  our 
titanium  bolometer  would  have  a  relatively  low  emissivity.  Because  the  bolometer  is 
typically  biased  with  low  voltages,  it  seems  unlikely  that  emissivity  would  play  a  role  at 
correspondingly  small  temperature  changes  (e.g.  less  than  100  K),  but  it  certainly 
warrants  analysis. 

The  question  is  whether  or  not  enough  power  is  radiated  away  fi*om  the  bolometer 
(Pradiated)  SO  as  to  make  a  significant  impact  on  the  measured  signal.  To  help  describe  the 
radiative  transfer  of  energy,  we  define  thermal  conductivity  due  to  radiation  Grad  as 
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p 

_  radiated 

rad- 


(2.21) 


with  the  understanding  that  it  is  not  an  actual  physical  parameter  but  rather  a  construct 
which  can  assist  us  in  describing  radiation  in  convenient  terms.  Substituting  equation 
(11.20)  for  the  radiant  power  into  equation  (2.21), 


G^AT  =  Ascr(T^-I^). 

Substituting  To+AT  for  T,  and  moving  a  dr  to  the  right  side, 

_Aea[lT,+AT]' -T,) 


(11.22) 


(11.23) 


Assuming  AT  is  small  relative  to  To* ,  all  but  the  first  order  terms  can  be  dropped  in  a 


Taylor  expansion  of  [2^  +  AT : 


[2;+Ar]*  =  r;  i+^ 

■^0  _ 

After  the  expansion  and  dropping  of  higher-order  terms, 

4Ar 


_-l4 


[r„+Ar] 


1  +  - 


This  approximation  can  be  re-inserted  into  equation  (11.23)  to  yield 


^rad  =■ 


Ascr 

( 

To 

1 

L  ^0  JJ 

AT 


(11.24) 


(11.25) 


(11.26) 


which  results  in 


G,^=^4AsaT,\  (11.27) 

This  result  surprised  the  author,  but  upon  reflection  it  certainly  makes  sense.  Every 

object  at  room  temperature  is  radiating  in  accordance  with  the  Stefan-Boltzman  Law;  a 

slight  change  in  temperature  will  not  have  a  major  effect  on  that  radiation.  The  fact  that 


*  At  low  bias  voltages,  typical  AT’s  are  around  30  °K,  compared  with  room  temperature  of  300  "K. 
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with  small  temperature  changes  Grad  is  constant  is  fortunate  for  those  who  wish  to 
establish  a  model  that  includes  it.  The  method  described  earlier  for  calculating  G  remains 
valid,  with  the  imderstanding  that  G  =  • 


6.  Figures  of  Merit 

The  figures  of  merit  for  a  general  thermal  detector  discussed  in  chapter  11. A.  1 
may  now  be  revised  so  as  to  describe  the  bolometer  in  a  Wheatstone  bridge  circuit.  We 
now  define  voltage  responsivity  as 


V 

R.=-^ 


B. 


(11.28) 


and  firom  equation  (II.2)  we  have 


AT  = 


incident 


p 

^  incident 


for  a  saturated  bolometer,  which  is  what  would  be  measured.  Substituting  equation  (II.2) 
into  equation  (II.  1 1),  we  have 


y'  _  '  b  ^  ^  incident 

4-G 


(11.29) 


and  equation  (11.28)  becomes 


Ry  =-^ — 

4-G 


F/ 


'W 


(11.30) 


Regarding  speed,  the  same  relationship  holds  true  for  the  time  constant  x,  which,  as  given 


by  equation  (II.4),  is 


7.  Measuring  Incident  Infrared  Radiation 

In  practice,  if  we  were  using  our  single  bolometer  Wheatstone  bridge  circuit  to 

make  measurements,  we  would  expose  it  to  the  IR  radiation  for  a  more  than  three  time 

constants  in  order  to  achieve  steady  state,  apply  a  bias  voltage  pulse  and  measure  the 

ouq)ut  voltage.  If  we  compared  this  output  voltage  data  with  that  of  an  experiment  in 
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which  we  just  measured  the  output  voltage  without  incident  IR,  we  would  see  something 
like  the  figure  given  below. 


Figure  II.  1 1  Exaggerated  illustration  of  the  output  voltage  to  due  Joule  heating  alone 
and  to  Joule  heating  combined  with  incident  IR. 

Figure  II.  11  illustrates  the  resultant  relative  output  voltages  and  is  only  provided  for 
pedagogical  use.  It  is  in  fact  a  gross  exaggeration;  the  difference  between  the  two  curves 
would  not  be  discemable  if  the  two  curves  could  be  drawn  to  scale.  Typical  values  for 
incident  IR  fluxes  are  in  the  nanowatt  (nW)  range,  compared  to  milliwatts  (mW)  due  to 
Joule  heating.  In  an  actual  circuit,  the  minuscule  IR  signal  would  have  to  be  extracted 
from  the  Joule  heating  response  prior  to  amplification.  If  it  were  not  removed, 
introduction  of  the  Joule  heating  signal  into  the  amplifier  would  cause  quick  amplifier 
saturation. 

One  way  to  compensate  is  to  electronically  subtract  a  calculated  value  of  the  Joule 
heating  voltage  from  the  composite  signal.  This  method  is  used  with  numerous  current 
devices,  but  it  relies  on  analytically  derived  values  to  be  subtracted  from  the  experimental 
composite  IR  signal  and  Joule  heating  signal.  If  the  signal  is  integrated  for  a  certain 
period  of  time  in  order  to  improve  the  S/N  ratio,  cancellation  becomes  even  more 
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difficult  since  the  self-heating  signal  increases  with  time  (see  equation  11.20)  until  the 
steady  state.  Another  way  would  be  to  physically  cancel  the  unwanted  voltage  by  using  a 
compensating  bolometer  on  the  opposite  side  of  the  Wheatstone  bridge,  as  illustrated  in 
the  following  schematic: 


iR 

??? 


Figure  II.  1 2  Compensating  bolometer  Wheatstone  bridge  circuit. 

In  this  arrangement  either  bolometer  Bi  or  bolometer  B2  is  the  “compensating” 
bolometer,  which  would  mean  that  it  should  be  optically  masked  from  the  incident  IR 
radiation,  but  allowed  to  imdergo  normal  Joule  heating.  With  the  objective  of  creating  a 
micro-machined,  integrated  circuit  on  a  silicon  substrate,  physically  shielding  the 
compensating  bolometer  from  the  IR  radiation  becomes  impractical.  However,  by 
altering  the  thermal  characteristics  of  the  compensating  bolometer,  it  can  be  made  to 
appear  virtually  translucent  to  the  IR  irradiation.  Specifically,  by  increasing  its  thermal 
conductivity,  i.e.  making  the  legs  which  attach  the  bolometer  face  to  the  heat  sink  larger, 
the  energy  which  is  absorbed  will  pass  through  the  bolometer’s  legs  to  the  heat  sink  at  a 
higher  rate. 
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8.  Compensating  Bolometer  Operations 


The  use  of  a  compensating  bolometer  with  higher  thermal  conductivity,  as 
opposed  to  an  identical  bolometer  which  is  physically  shielded  from  incident  IR  flux, 
introduces  a  host  of  dynamic  changes  to  the  operation  of  the  Wheatstone  bridge. 
Changes  in  signal  strength  as  well  as  thermal  time  constants  accompany  any  change  in 
the  thermal  conductivity  of  the  reference  bolometer.  Figure  11.13  depicts  the  temperature 
changes  experienced  by  a  sensor  bolometer  Bi  (upper  curve)  and  a  reference  bolometer 
B2  (lower  curve)  as  they  are  both  exposed  to  incident  IR  flux  and  then  to  a  bias  voltage 
step  function. 


Figure  11.13  Temperature  responses  of  a  sensor  bolometer  and  reference  bolometer. 
The  upper  curve  represents  Bi  and  the  lower  curve  represents  B2,  where  G2=2Gi. 


In  the  interest  of  simplicity,  incident  IR  power  has  been  made  equal  to  the  power 
dissipated  due  to  bias  voltage.  The  thermal  conductivity  of  the  reference  bolometer  has 


26 


been  arbitrarily  selected  to  be  twice  that  of  the  sensor  bolometer,  in  order  to  demonstrate 
how  thennal  conductivity  affects  both  the  maximum  signal  strength  and  the  speed  with 
which  the  bolometer  heats  up,  in  accordance  with  equations  (II.2)  and  (11.4),  respectively: 


II 

<1 

(11.2) 

> 

II 

1-e‘^'j  where  t  = 

(11.4) 

The  signal  we  wish  to  measure  is  the  maximum  temperature  of  the  sensor 
bolometer  due  to  incident  IR  flux,  i.e.,  the  upper  curve  before  bias  voltage  is  applied. 
With  the  reference  bolometer  in  place,  the  actual  IR  signal  available  is  the  difference 
between  the  two  curves  before  the  bias  voltage  (AT).  From  this  perspective,  the  reference 
bolometer  has  degraded  the  signal  (by  half).  In  practice,  if  we  use  a  reference  bolometer 
with  a  thermal  conductivity  one-tenth  that  of  the  sensor  bolometer,  we  will  experience 
only  a  10%  signal  loss.  Hence,  the  first  of  two  design  constraints  is  evident:  G2  should  be 
as  large  as  possible  to  make  B2  optically  transparent  to  IR  radiation. 

The  benefit  of  the  reference  bolometer  becomes  manifest  when  the  bias  pulse  is 
applied  and  the  two  temperature  levels  rise  exponentially.  Considering  the  more  realistic 
case  in  which  Joule  heating  power  is  orders  of  magnitude  higher  than  IR  power,  one 
might  be  tempted  to  wonder  how  an  accurate  signal  can  be  extracted  at  all  with  the  use  of 
a  reference  bolometer.  Analytic  results  in  the  case  where  G2=2Gi  and  dissipated  power 
equals  1000  times  incident  IR  power  are  given  in  the  following  figure,  generated  by 
MATLAB. 
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Figxire  11. 14  Scale  view  of  response  in  which  applied  power  =  lOOOxIR  power  and 
G2=2xGi.  The  right  hand  figure  illustrates  the  initial  equivalent  linear  growth  for  both 

curves. 

The  original  signal  is  completely  buried  shortly  after  bias  voltage  is  applied  for  the  read¬ 
out.  Fortunately,  while  the  temperatures  of  Bi  and  B2  grow  with  different  exponential 
rates,  there  is  a  region  in  which  their  growth  rates  may  be  approximated  as  linear  and 
independent  of  thermal  conductance,  which  is  illustrated  in  the  right  hand  graph  of  Figure 
11.14.  Rearranging  the  heat  transfer  equation, 

=  .  (11.31) 

H  T 

If  dt«.T ,  then  equation  (II.3 1)  can  be  approximated  as 

d(M)=^^dt.  (11.32) 

Therefore,  in  the  region  dt«.r,  where  dt  is  the  bias  voltage  pulse  width,  we  should  be 
able  to  integrate  the  signal  to  measure  AT  quite  accurately,  as  long  as  the  heat 
capacitances  of  the  two  bolometers  are  the  same.  Here  we  see  om  second  design 
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constraint:  for  x  to  be  large  compared  to  dt,  G2  must  be  small,  so  that  it  may  heat  up  at 
the  same  rate  as  Bi  for  at  least  long  enough  to  obtain  a  read-out  from  integration. 

These  two  competing  design  constraints  were  the  motivation  for  a  computer 
model  that  could  be  used  to  theoretically  optimize  various  Wheatstone  bridge 
configurations.  With  this  in  mind,  Section  III  deals  with  developing  and  verifying  an 
acceptable  model  using  an  industry-standard  electrical  engineering  tool  called  PSPICE. 
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III.  MODELING  THE  BOLOMETER  WITH  PSPICE 


A.  DISCUSSION  OF  PSPICE 

PSPICE  is  a  computer  simulation  program  that  models  the  behavior  of  an 
electrical  circuit  containing  analog  devices.  It  is  essentially  a  computer-based  breadboard 
which  allows  one  to  predict  AC  and  DC  steady  state  and  transient  analyses.  By  using  it 
to  model  the  circuits  described  in  Section  II  and  then  verifying  the  model  against 
experimental  data,  we  may  confidently  conduct  parametric  studies  of  the  bolometer. 

1.  Bolometer  Heat  Transfer  Equation  versus  the  Parallel  RC  Circuit 

As  discussed  in  Section  I,  the  heat  transfer  equation  is  of  the  same  form  as  the 
differential  equation  for  the  parallel  RC-circuit  under  constant  current  charging  (Figure 
II.3.)  The  following  are  the  two  analogous  equations: 

P(r)  =  g»Ar  +  (P(/)  constant) 

dt  ^  (III.1) 

/(r)  =  i(F(t)-Fo)  +  C^^  (/(?)  constant). 


Their  analogous  parameters  are  as  follows,  with  their  SI  units  in  square  brackets: 


T{t)  [Ji:]  ^  v(t)  [V],  (ni-2) 

p{t)  [w]  -*  I{t)  [.4]. 


Using  this  analogy,  we  can  design  a  circuit  using  PSPICE  to  mimic  the  temperature- 

dependency  of  the  bolometer’s  resistance.  We  can  then  insert  this  virtual  circuit  into 
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another  circuit  representing  the  actual  Wheatstone  bridge  with  a  bias  voltage  applied  to  it. 
Special  consideration  must  be  given  not  to  confuse  the  actual  electric  circuit  values  with 
the  virtual  “temperature  circuit”  values.  The  bolometer  is  modeled  as  an  RC  circuit  with 
the  resistor  representing  the  thermal  conductance  and  the  capacitor  representing  heat 
capacity,  as  shown  below. 


p  —  bolometer^ 

-  /R{T) 


bolometer 


Cy=H, 


Figure  III.  1  Virtual  circuit  representing  a  bolometer. 


Subscripts  “R”  and  “V”  stand  for  real  and  virtual  quantities,  respectively.  R(T)boiometer  is 
given  by  R(T)  =  Rofl +a-AT),  and  Vboiometer  is  the  value  of  the  voltage  drop  across  the 
bolometer  in  the  Wheatstone  bridge.  The  value  of  the  current  source  in  the  virtual  circuit 
corresponds  to  the  power  dissipated  in  the  bolometer.  In  figure  III.l,  a  voltage  source 
maintains  a  virtual  offset  voltage  of  300  V,  placing  the  sensor’s  initial  temperature  at 
room  temperature  conditions  of  300  K.  What  remains  is  to  connect  the  bolometer  to  a 
Wheatstone  bridge  and  apply  a  bias  voltage,  all  of  which  we  can  do  using  PSPICE. 
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B,  THE  SINGLE  BOLOMETER  WHEATSTONE  BRIDGE 

1.  PSPICE  Model  of  Single  Bolometer  Wheatstone  Bridge 

The  point  of  interface  for  the  circuit  described  above  is  the  bolometer  itself, 
represented  by  an  analogue  behavior  modeling  (ABM)  box,  which  allows  the  user  to 
control  the  output  through  the  use  of  the  algorithms  he  chooses.  As  there  is  no 
temperature-dependent  resistor  in  the  ABM  options  provided  by  PSPICE,  we  must 
replace  the  function  of  the  resistor  with  an  available  option:  a  variable  current  source, 
which  we  can  program  to  set  the  current  through  to  be  the  same  as  that  through  the 
bolometer  it  replaces. 


Figure  III.2  PSPICE  model  of  Wheatstone  bridge  with  a  single  bolometer. 


Figure  111.2  represents  the  circuit  as  seen  in  the  PSPICE  format,  with  the  exception  being 
that  the  ABM  formula  accept  only  numeric  values,  as  opposed  to  the  formulas  shown  in 
the  figure,  which  include  a  and  Ro  to  help  clarify  the  figure.  The  ABM  box  receives 

temperature  inputs  — from  the  virtual  circuit  to  insert  the  correct 


33 


cxxrrent  into  the  Wheatstone  bridge,  acting  as  if  there  were  a  temperature-dependent 
resistor  there.  The  current  generated  is  given  by 

T ( t\  —  ^bolometer /  ,  , 

/Riry 

where  R(T)  is 

i?(7’)  =  i?o(l  +  a.Ar). 

The  virtual  circuit  is  able  to  calculate  the  temperature  difference  AT  by  inserting  the 
voltage  drop  across  the  bolometer  in  the  actual  circuit  into  an  ABM  box  (far  left)  to 
generate  a  virtual  power,  whose  RC-circuit  analogue,  according  to  equation  (III.2),  is 
current. 

At  this  point  the  circuit  has  been  designed  to  represent  just  the  Joule  heating  effect 

without  incident  radiation.  To  include  irradiance,  a  current  source  (analogous  to  power) 

would  be  placed  in  parallel  with  the  resistor  and  capacitor  in  the  bolometer  circuit. 

2.  Verification  of  the  PSPICE  Model  and  Early  Comparison  with  the 
Analytic  Model 

After  creating  the  PSPICE  model  shown  in  figure  III.2,  the  first  step  to 
verification  is  to  observe  the  correlation  of  simulated  data  to  the  analytical  predictions  of 
Section  II.  One  of  the  first  experiments  run  was  to  verify  equation  (11.13),  which  states 
that  the  dissipated  power  in  the  bolometer  is  constant,  given  a  steady  bias  voltage  (F*): 

(n.i3) 


(III.3) 


(111.4) 
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Recalling  that  a  quadratic  term  of  temperature  (R^a^AT^)  was  dropped  from  the 
denominator,  what  this  simulation  really  seeks  to  determine  is  for  what  region  equation 
(11.13)  is  valid.  Using  a  bolometer  with  the  same  values  as  before  (G  =  3.2x10"’  IF/ 

H  =  J ! K ,  R  =  \6kO.,  and  a=0.0027),  the  results  for  applied  bias  voltages 
ranging  from  1.0  to  10.0  volts  are  given  below. 


Figure  III.3  Joule  heating  power  dissipated  for  varying  bias  voltages. 


Also  shown  are  the  corresponding  analytical  values  for  the  top  three  bias  voltages  (8.0, 
9.0, 10.0  V)  as  given  by  equation  (11.13).  What  is  quite  apparent  is  the  departure  from  the 
expected  linear  relationship.  The  nonlinear  effects  of  the  disregarded  quadratic  and 
higher  temperature  terms  are  in  fact  important.  So  what  is  the  acceptable  region  of 
linearity? 
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We  can  determine  the  condition  for  acceptably  disregarding  the  quadratic  term  in 


the  equation 


P  =  - 


R^{\  +  a^T)  ^^2 
4J^+4Ii^aAT+Ii^a^Ar^  * 


(IIL5) 


by  setting  up  the  condition  to  be 


4Ji^aAT»Ji^a^AT^ 


which  then  leads  to 


-»Ar.  (III.6) 

a 

In  the  case  of  a=0.0027  for  Ti  bolometers,  the  acceptable  region  must  be  AT  1480  K. 
The  exact  region  obviously  depends  upon  how  accurate  a  prediction  one  would  like  to 
make.  The  deterioration  in  linearity  can  clearly  be  seen  in  the  following  set  of  graphs. 


Time  (ms)  Time  (ms) 
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Power  Dissipated  (mW)  Power  Dissipated  (mW)  Power  Dissipated  (mW) 


Power/Temperature  Relationship  for  bias  Voltage  =  5.0  V 


Power/Temperature  Relationship  for  bias  Voltage  =  6.0  V 


Power/Temperature  Relationship  for  bias  Voltage  =  7.0  V 


Figure  III.4  Effects  of  increasing  temperature  on  dissipated  power. 
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One  can  actually  begin  to  see  degradation  on  the  second  graph  (Vb=2.0  V),  and  by  the 
fourth  graph  (Vt,=4.0  V)  there  is  significant  deviation  from  the  expected  linear  behavior. 
At  Fi=10.0  V,  the  steady  state  power  is  slightly  over  half  of  the  value  predicted  by  the 
analytic  model.  By  this  point  the  temperature  change  is  at  approximately  1600  K.  Hence 
we  have  clear  verification  that,  in  order  for  our  model  to  be  accurate,  the  bolometer  must 
be  operated  in  the  region  well  below  1480  K,  as  equation  (111.6)  states. 

3.  Verification  of  the  Model  by  Comparison  with  Physical  Data 

Having  established  the  acceptable  temperature  regime,  we  can  also  verify  the 
model  by  comparison  with  actual  physical  data  (Ramakrishna,  et  al,  1999).  Two 
bolometers  with  distinct  thermal  properties  were  used  in  a  preliminary  experiment  to 
demonstrate  the  concept  of  a  compensating  bolometer  Wheatstone  bridge.  The  thermal 
characteristics  are  given  in  the  following  table. 


Bolomter  Bi 

Bolomter  B2 

Resistance  (at  294  K) 

26  kQ 

24  kQ 

a  (at  294  K) 

0.0027/K 

0.0029/K 

Thermal  Conductance  (G) 

3.61x10’’ W/K 

2.37x10’^  W/K 

Thermal  Mass  (H) 

4.0x10’*  J/K 

4.74x10’*  J/K 

Thermal  Time  Constant  (x) 

111  ms 

15  ms 

Responsivity  (at  3  V) 

6000  V/W 

770  VAV 

T able  HI.  1  Bolometer  characteristics 


Below  is  a  graph  showing  measured  output  voltage  versus  time  for  the  two 
bolometers.  Note  that  the  upper  and  lower  x-axes  separately  represent  the  two 
bolometers. 
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Time  axis  for  Bj  (ms) 

0  20  40  60  80  100 


Figure  III.  5  Experimental  output  voltages  for  the  two  bolometers  as  a  function  of  time 

at  294  K  (Ramakrishna,  et  al). 


Running  a  PSICE  simulation  with  the  same  parameters,  the  following  data  is  generated: 


Time  axis  for  (ms) 

0  20  40  60  80  100 


Figure  111.6  Simulated  output  voltage  for  bolometer  Bi  and  bolometer  B2  at  294  °K. 


The  agreement  between  Ramakrishna  and  Karunasiri’s  physical  experiment  and  the 
simulated  experiment  of  figure  111.6  falls  within  1%.  The  data  for  figure  111.6  is  collected 
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in  PSPICE  by  placing  “voltage  differential  probes”  on  either  side  of  the  Wheatstone 
bridge,  in  each  case  between  the  two  resistors,  as  shown  in  figure  III.2.  We  can  simulate 
monitoring  the  temperature  by  placing  the  voltage  differential  probes  on  the  upper  and 
lower  branches  of  the  bolometer’s  virtual  temperature  circuit  in  figure  III.2.  Doing  so 
yields  the  following  data. 


Time  axis  for  (ms) 

0  20  40  60  80  '100 


Figure  III.7  PSPICE  simulated  temperature  change  for  bolometer  Bi 

While  it  is  not  possible  to  measure  the  steady-state  temperature  change  in  the  physical 
experiment,  it  can  at  least  be  calculated  anal)dically  to  determine  how  closely  the  model 
agrees  with  analytic  predictions.  For  bolometer  Bi,  with  Vout=  0.0091  V,  Vb=  0.8  V  and 
a  =  0.0027  ,  equation  (H.l  1)  yields  16.85  K.  The  PSPICE  simulation  yielded  a 

maximum  temperature  change  of  17.02  K,  which  indicates  agreement  within  1%.  For 
bolometer  B2,  with  Voui=  0.007  V,  Vb=L3  V  and  a  =  0.0029,  equation  (II.  1 1)  yields  AT= 


40 


7.427  K.  In  this  case  the  PSPICE  simulation  generated  a  steady-state  temperature  change 
of  7.416  K,  well  under  a  1%  discrepancy.  It  should  be  noted  that  these  slight  errors  may 
be  due  to  estimating  the  saturation  voltages  from  figure  III.6,  as  actual  numerical  data 
was  not  available  from  the  physical  experiment. 

Before  completing  the  compensating  bolometer  Wheatstone  bridge,  we  can 
compare  modeled  data  with  the  physical  data  on  how  the  individual  bolometers  respond 
to  Joule  heating  and  incident  IR  in  the  dt«.z  region,  where  dt  is  the  bias  pulse  width. 


Figure  II1.8  Physical  data  for  bolometer  Bi  (left)  and  bolometer  B2  (right). 

Once  again,  using  bolometers  with  identical  thermal  characteristics,  the  following  output 
voltage  data  is  generated  by  PSPICE.  (next  page) 


41 


Figure  III.9  PSPICE  simulation  of  bolometer  response  to  incident  IR  in  addition  to 
Joule  heating.  Bolometer  Bi  is  represented  by  the  left  graph,  bolometer  B2  by  the  right 

graph. 


In  both  figures  III.8  and  III.9,  we  can  see  the  IR  signal  riding  on  top  of  the  Joule 
heating  signal.  Additionally,  the  IR  signal  is  significantly  stronger  on  bolometeri  than  on 
bolometera,  which  was  the  goal  of  attempting  to  make  bolometer2  translucent  to  IR  by 
increasing  its  thermal  conductance.  Nevertheless,  because  we  are  operating  in  the 
dt<^z  region,  the  output  voltages  due  to  Joule  heating  for  both  bolometers  are  roughly 
equivalent,  following  the  theory  leading  to  equation  (11.32).  The  PSPICE  model  appears 
quite  accurate  even  at  small  transient  times. 


4.  Modeling  the  Compensating  Bolometer  Wheatstone  Bridge  and  Final 
Verification. 

With  the  pieces  in  place  and  fimctioning  properly,  it  remains  to  assemble  the 
Wheatstone  bridge  with  a  compensating  bolometer  in  place  in  order  to  cancel  out  the 
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majority  of  the  Joule  heating  signal,  leaving  the  IR  signal  and  a  roughly  equivalent  Joule 
heating  signal  as  the  input  to  the  read-out  circuitry,  which  would  ostensibly  begin  with  an 
amplifier  (the  very  reason  why  it  is  advantageous  to  cancel  out  the  enormous  Joule 
heating  signal).  The  following  physical  data  show  good  cancellation  of  the  Joule  heating 
signal. 


Time  (ms) 

Figure  III.  1 0  IR  and  Joule  heating  response  of  the  bridge  consisting  of  two  bolometers 
contrasted  with  response  of  the  bolometer  Bi  at  bias  voltage  of  3V.  The  two  bolometers 
were  exposed  to  the  same  amoimt  of  infrared  from  a  blackbody  source. 

The  actual  numerical  physical  data  was  not  available,  nor  were  the  specifics  of  the 
experiment  setup.  For  the  PSPICE  simulation,  the  Wheatstone  bridge  was  set  up  with  Bi 
on  one  side  with  an  equivalent  resistor  and  B2  on  the  other  side  with  an  equivalent 
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resistor.  (This  setup  is  shown  in  figure  III.  12.)  The  simulation  yielded  the  following 
data. 
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Figure  III.  1 1  PSPICE  simulation  of  the  IR  and  Joule  heating  response  of  the  bridge 
consisting  of  two  bolometers  contrasted  with  response  of  the  bolometer  Bi  at  bias  voltage 
of  3  V.  The  two  bolometers  were  exposed  to  the  same  amount  of  simulated  infrared  flux 

from  a  blackbody  source. 


The  correlation  between  the  physical  data  and  the  simulation  data  is  favorable, 
once  again  under  1%  difference.  Nevertheless,  while  the  comparison  of  the  physical 
cancellation  with  the  simulated  cancellation  are  inconclusive,  the  individual  bolometers 
matched  very  well,  and  it  is  assumed  that  combining  the  two  onto  one  compensated 
Wheatstone  bridge  should  not  alter  their  individual  performance,  as  each  is  subject  to  the 
same  voltage.  Therefore,  on  the  basis  of  the  excellent  correlation  between  reality  and 
simulation  as  demonstrated  by  figures  III.  8  and  III.9,  the  model  is  deemed  a  success. 
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Figure  III.  1 2  Compensated  bolometer  Wheatstone  bridge 

[Note:  to  include  variable  thermal  radiation  for  cases  when  the  condition  AT  <§:  does  not  apply,  ABM 
boxes  should  be  generated  in  parallel  in  both  bolometer  circuits  and  programmed  with  the  Stefan-Boltzman 
Law,  equation  (11.20).] 
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c. 


PARAMETRIC  STUDY  USING  THE  WHEATSTONE  BRIDGE  MODEL 


The  virtue  of  the  PSPICE  model  can  now  be  realized  in  that  it  is  readily  possible 
to  conduct  research  and  design  studies  of  microbolometers  without  actually 
manufacturing  them,  saving  time  and  money.  Given  a  certain  application,  the  designer 
may  now  quickly  produce  the  necessary  thermal  characteristics  required  to  design  a 
useable  device.  It  should  be  cautioned,  however,  that  the  assumptions  and 
approximations  made  in  designing  the  model  must  still  hold  true  in  order  to  create  an 
accurate  design.  What  follows  is  an  example  of  how  bolometer  applications  can  be 
improved  by  employing  the  compensation  technique,  as  illustrated  through  the  use  of  the 
PSPICE  model.  The  parameter  to  be  used  for  comparison  will  be  integration  time,  which 
therefore  warrants  a  brief  discussion. 

1.  Integration  Time  versus  Noise  Bandwidth 

A  cursory  look  at  the  forms  of  noise  affecting  bolometers  reveals  that  they  are 
dependent  on  frequency  bandwidth  Af: 


=  AkTR  ■  Af  Johnson  noise 


where  R  is  resistance,  k  is  Boltzman’s  constant,  T  is  temperature,  and  5  is  a  constant  for 


noise.  (Shot  noise  is  deemed  negligible  here).  In  our  application,  Af  is  determined 


by  the  pulse  width  of  the  sampling  pulse,  where  pulse  width  directly  correlates  to 
integration  time  tint.  Recalling  from  Fourier  analysis,  a  narrower  square  pulse  requires  the 
addition  of  higher  frequencies  for  its  construction  than  a  broader  square  pulse.  Generally, 
the  bandwidth  is  given  by 
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Af — (ni.7) 

Hence,  the  longer  one  can  integrate  a  signal,  the  less  noise  will  remain.  In  the 
comparison  which  follows,  the  decision  is  arbitrarily  made  that  the  signal  will  be 
integrated  for  a  time  period  sufficient  for  the  voltage  response  due  to  Joule  heating  to 
equal  the  voltage  response  due  to  an  arbitrarily  chosen  IR  radiation  flux.  This  concept 
reflects  the  feasibility  that  the  subsequent  electronic  subtraction  of  a  Joule  heating  voltage 
on  the  same  order  of  magnitude  as  the  signal  voltage  reduces  the  probability  of  distorting 

the  remaining  signal  voltage. 

2.  Application:  Thermal  Imaging  Device 

The  necessity  and  usefulness  of  lightweight,  portable  thermal  imaging  devices  are 
obvious  to  military,  law  enforcement  and  fire  fighting  personnel.  The  ability  to  image 
targets,  suspects,  and  fires  or  hot  spots  in  completely  dark  and/or  smoky  scenarios  is 
critical  to  the  success  of  many  kinds  of  operations  carried  out  by  such  personnel.  As 
explained  in  the  introduction,  there  are  currently  numerous  devices  available,  but 
typically  the  operator  must  accept  a  trade-off  in  terms  of  sensitivity  versus  weight, 
because  the  more  sensitive  devices  typically  need  to  be  be  cooled,  which  requires  more 
power,  weight  and  volrnne.  The  microbolometer  fills  this  niche  of  high  sensitivity  and 
low  weight/space  at  the  expense  of  speed.  Bolometers  are  relatively  slow  devices 
compared  with  many  others,  but  this  lack  of  speed  really  only  precludes  them  firom  high¬ 
speed  applications,  such  as  IR  warhead  seekers. 

Raytheon  Systems  Company  has  developed  a  prototype  infirared  imaging  rifle- 
sight  using  an  imcooled  microbolometer  focal  plane  array  (FPA)  (Radford,  et  al,  1998). 
(Implemenatation  of  the  bolometer  into  a  FPA  is  introduced  in  Appendix  A.)  This 
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system  is  typical  for  the  kind  of  application  needed  in  the  Department  of  Defense. 
Because  it  is  an  integrated,  micromachined  chip  technology,  it  is  relatively  low  in  cost 
compared  with  other  microbolometer  applications  that  require  complicated 
manufacturing  procedures  to  hide  the  reference  bolometer  from  IR  radiation,  for 
example.  In  this  case,  the  elemental  structure  is  a  one-bolometer  Wheatstone  bridge,  and 
the  Joule  heating  voltage  is  electronically  subtracted  using  a  dark  current  subtraction 
method.  The  thermal  characteristics  of  the  bolometer  are  given  below. 


Performance  Parameter 

Capability 

Pixel  Dimensions 

50  (pm)^ 

Thermal  Conductance  (G) 

<5xl0'*W/K 

Thermal  Capacity  (H) 

<1x10-^  J/K 

Thermal  Time  Constant  (x) 

20  msec 

a 

>.02 

Resistance  {Rq) 

50  kQ 

Table  III.2  Ra5dheon  IR  imaging  rifle-sight  bolometer  characteristics. 

Given  such  a  bolometer,  we  would  like  to  determine  the  advantage  gained  in  using  a 

compensating  bolometer  in  the  Wheatstone  bridge  configuration.  From  equation  (11.15) 

we  see  that  the  saturation  voltage  from  a  one-bolometer  Wheatstone  bridge  is  given  by 


V.-aP 

4  G' 


(III.8) 


However,  we  must  consider  the  saturation  voltage  for  a  Wheatstone  bridge  with  an 

additional  compensating  bolometer,  whose  thermal  conductance  is  different  from  the 

sensor  bolometer  by  definition.  This  calculation  would  be  analytically  impossible  were  it 

not  for  the  principle  of  superposition  (Boas,  1983),  which  allows  us  to  analyze  the 

Wheatstone  bridge  twice,  once  each  with  a  different  bolometer,  and  then  add  the  results. 

The  principle  of  superposition  is  valid  for  any  linear  network  containing  resistors  and 
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voltage  sources.  Using  this  method  in  the  situation  where  the  compensating  bolometer  is 
identical  to  the  sensor  bolometer  with  the  exception  of  thermal  conductance,  we  retrace 
the  steps  which  lead  to  equation  (11.15)  and  arrive  at 


F,-a-P  1 _ 

4  [g,  G2 


(111.9) 


A  convenient  figure  of  merit  for  the  detector  would  be  its  signal  response  at  4  nW  of 
incident  IR  flux.  For  the  Raytheon  bolometer  circuit  the  signal  would  be  1.2  mV  at  Vb=3 


V.  For  a  reference  bolometer  with  thermal  conductance  G2  one  order  of  magnitude 


greater  (5x10'’  W/K),  the  signal  would  be  1.08  mV.  There  has  been  a  10%  reduction  in 
signal  strength,  and  while  the  signal  is  certainly  strong  enough  to  manipulate,  any 
reduction  in  signal  strength  is  not  desirable.  However,  the  following  graph  illustrates  the 
gain  in  integration  time,  assuming  the  signal  will  be  integrated  until  the  Joule  heating 
matches  the  signal  strength. 


Time  (jis)  Time  (jis) 


Figure  III.  1 3  Comparison  of  integration  times  with  and  without  a  compensating 

bolometer. 
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Clearly  the  benefit  in  integration  time  far  outweighs  this  loss  of  signal  strength,  when  one 
considers  the  effect  on  noise  that  an  integration  time  of  approximately  two  microseconds 
versus  approximately  ninety  microseconds  has.  Following  equation  (111.7),  the  one 
bolometer  circuit  has  a  bandwidth  of  5x10^  Hz  while  the  compensated  circuit  has  a 
bandwidth  of  1.1  x  10"^  Hz.  It  follows  that  the  noise  in  the  compensated  bolometer  is  one- 
and-a-half  orders  of  magnitude  less  than  in  the  uncompensated  circuit. 

For  imaging  devices,  this  kind  of  noise  reduction  greatly  enhances  image  quality. 
Video  imagery  is  typically  around  30  frames  per  second  (fps),  and  even  with  the 
extended  integration  time  we  are  nowhere  near  this  limit.  Even  with  the  10%  loss  of 
signal,  implementation  of  a  compensating  bolometer  in  the  IR  imaging  rifle-sight  would 
improve  it  dramatically. 

3.  Further  Possible  Benefit  -  Cancellation  by  Integration 

The  benefit  of  having  developed  a  working  model  of  the  compensated 
Wheatstone  bridge  becomes  manifest  when  one  simply  begins  to  tinker  with  the  design  of 
the  compensating  bolometer.  Instead  of  spending  valuable  resources  in  the  design, 
manufacture  and  testing  of  physical  bolometers,  concepts  can  be  modeled  with  ease. 
Such  simulations  may  hold  the  key  to  the  actual  elimination  of  the  Joule  heating  output 
voltage.  In  adjusting  the  thermal  capacities  and  conductivities  somewhat  haphazardly,  it 
was  found  that  a  Joule  heating  curve  could  be  generated  which  first  went  negative,  and 
then  curved  back  positive  before  it  exploded  into  the  normal  exponential  growth  cxirve. 
This  was  intriguing  because  it  appeared  to  offer  a  means  of  canceling  out  the  Joule 
heating  voltage;  simply  integrating  for  a  period  long  enough  to  cancel  the  negative 
voltage  area  with  an  equivalent  positive  voltage  area  would  eliminate  the  Joule  heating 
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voltage.  Unfortunately,  it  was  noted  that  adjusting  the  thermal  mass  of  the  reference 
bolometer  would  be  very  impractical  in  the  micro-machining  process,  and  therefore 
prohibitively  expensive. 


Fortunately,  from  the  analytic  description  of  the  bolometer’s  early  time  response 
another  option  for  achieving  the  same  behavior  appears: 


V 

16R. 


H 


■t. 


(11.19) 


Equation  (11.19)  indicates  that  output  voltage  is  equally  dependent  on  the  balancing 
resistor  values  as  the  thermal  mass.  This  leads  to  the  idea  that  the  same  advantageous 
output  voltage  curve  can  be  generated  by  manipulating  the  value  of  Ro  on  one  side  of  the 
Wheatstone  bridge.  The  figure  below  illustrates  the  idea. 


Figure  III.  1 4  Unbalanced  Wheatstone  bridge 


The  design  shown  in  figure  111.14  was  implemented  in  PSPICE,  with  the  same 
Raytheon  bolometer  as  described  by  table  111.2.  This  time,  instead  of  precisely  balancing 
R]  and  R2  to  equal  each  other,  the  value  of  R2  was  adjusted  downward,  resulting  in  the 
following  set  of  ouQ)ut  curves. 
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20 


Figure  III.  1 5  PSPICE  simulated  output  voltage  versus  time  for  integration  by 

cancellation  design. 

As  hoped  for,  adjusting  R2  lower  than  Rj  caused  the  curve  to  dip  negative  before 
launching  exponentially  in  the  positive  direction.  To  integrate  the  Joule  heating  voltage, 
an  integrator  circuit,  such  as  a  Miller  integrator,  would  have  to  be  integrated  into  the  pre¬ 
amplifier  circuitry.  It  would  therefore  be  important  to  keep  the  negative  voltage  dip 
roughly  equivalent  to  the  signal  voltage,  in  order  to  keep  the  integrator’s  capacitor  small 
enough  for  an  appreciable  output  signal.  R2  meets  this  requirement,  while  providing  for 
an  integration  time  of  approximately  450  microseconds.  The  Joule  heating  voltage  could 
be  entirely  cancelled  out,  leaving  behind  only  the  signal  voltage.  Not  only  is  the  sensor 
greatly  enhanced  by  removal  of  the  Joule  heating  voltage,  but  the  integration  time  has 
been  extended,  further  improving  the  noise  reduction  by  integration,  and  hence  the 
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signal-to-noise  ration  (SNR).  For  example,  the  predominate  source  of  noise,  Johnson 
noise,  may  be  calculated  as  follows. 


Using  the  principle  of  superposition,  we  can  derive  the  noise  generated  by  each  of 
the  four  resistors  individually,  and  then  add  the  noise  figures  in  quadrature,  as  should  be 
done  with  uncorrelated  noise  sources.  Assuming  all  noiseless  resistors  with  one  noise 
voltage  source  for  the  lower  right  resistor,  the  circuit  appears  below. 

Vr 

Dias 


Figure  III.  1 6  Wheatstone  bridge  with  voltage  noise  source  for  lower  right  resistor 


The  development  is  similar  to  that  for  equation  (II.  1 1),  but  this  time  V2  is  different,  with 

V  -V 

(III.  10) 

instead  of  Fj  =  ,  where  V„  is  the  noise  voltage  due  to  one  resistor.  It  then  follows 


that 


AK, 


V 

AF  .  =-^. 

noise  ^ 


(III.  11) 
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For  the  total  noise  voltage,  the  individual  noise  voltages  from  each  resistor  are  added  in 
quadrature.  In  this  case,  the  four  individual  noise  voltages  are  equivalent,  so 


or 

In  dealing  with  Johnson  noise. 


(ni.l2) 


Vl,oJSE=^l^kTR^ .  (111.13) 

We  may  now  evaluate  the  benefit  of  the  longer  integration  time  provided  by  the 

cancellation  by  integration  technique.  From  figure  III.9  we  conservatively  estimate  an 

integration  time  of  t  =  2  ps,  which  also  leaves  behind  0.6  mV  of  Joule  heating  voltage. 


Using  this  integration  time,  R=  50  kQ  and  r=294  K,  equation  (III.  13)  yields 

^NoisE  ~  14.25  fjV .  (III.  14) 

Comparing  this  with  cancellation  by  integration  using  Rz  =  47kQ,  we  can  estimate  x  = 


425  ps  from  figure  II.l  1.  This  time  equation  (III.  13)  yields 


^NOISE  •  (III.  15) 

In  essence,  these  figures  represent  the  minimum  detectable  signal  threshold,  because  the 

sampling  time  is  still  short  enough  that  VC  noise  is  negligible  compared  to  Johnson 


noise.  In  addition,  the  normally  remaining  Joule  heating  voltage  can  be  completely 
eliminated,  fiirther  reducing  the  inaccuracies  that  accompany  electronic  subtraction  of  the 
Joule  heating  voltage. 


This  appears  to  be  a  very  promising  technique,  because,  while  altering  an 
individual  bolometer’s  thermal  mass  (H)  is  very  difficult  (requiring  a  three-dimensional 
architectural  change),  altering  its  thermal  conductance  (G)  and  baseline  resistance  is 
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simply  a  matter  of  changing  the  two-dimensional  architecture  of  the  reference  bolometer 
in  the  micro-machining  process. 
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IV.  CONCLUSION 

By  developing  a  working  model  on  PSPICE,  we  have  demonstrated  the  ability  to 
explore  more  deeply  the  complex  interaction  between  the  thermal  characteristics  of  the 
paired  bolometers  in  the  compensating  Wheatstone  bridge.  Understanding  the  complex 
thermal  interactions  is  only  possible  because  the  computer  model  numerically  accounts 
for  the  higher-order  terms  which  are  neglected  by  necessity  in  analytic  treatments.  In 
practical  terms,  the  model  allows  the  researcher  to  tinker  and  to  explore  these  interactions 
without  incurring  enormous  costs  or  wasting  valuable  research  time. 

One  example  of  the  benefit  of  this  kind  of  free  tinkering  was  arrived  at  very 
quickly,  with  the  discovery  of  cancellation  by  integration  while  using  an  unbalanced 
bridge.  Previous  thought  had  been  entirely  directed  at  the  balanced  bridge.  Through  the 
cancellation  by  integration  technique  developed  in  section  III.C.3,  we  were  able  to  realize 
a  dramatic  increase  in  integration  time,  which  equated  to  a  minimum  detectable  signal 
more  than  an  order  of  magnitude  smaller  than  without  the  technique. 

A  great  deal  of  information  awaits  those  who  wish  to  more  thoroughly  explore  the 
Wheatstone  bridge  arrangement  using  this  model.  With  five  independent  variables  all 
affecting  the  output  voltage,  there  remain  many  imusual  configurations  worthy  of  study. 
One  obviously  worthwhile  pursuit  would  be  to  study  the  sensitivity  to  manufacturing 
tolerances,  in  order  to  determine  how  tightly  a  particular  design  must  be  manufactured  to 
deliver  a  tolerable  performance.  Another  beneficial  pursuit  would  be  that  of 
understanding  the  physical  limits  of  a  microbolometer  based  on  the  architectural  design 
limits  of  the  micromachining  process  (i.e.,  how  narrow  a  bridge  can  be  made  and  what 
this  means  in  terms  of  thermal  conductance  limits). 
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APPENDIX  A:  THE  BOLOMETER  ARRAY 


A  bolometer  array  can  be  designed  with  the  reference  bolometer  outside  of  the 
pixel  array.  Using  vertical  and  horizontal  shift  registers,  each  individual  bolometer  can 
be  read  sequentially  and  incorporated  into  the  Wheatstone  bridge  structure  in  the  lower 
right  comer. 


Figure  A.  1  Microbolometer  Focal  Plane  Array  (FPA)  architecture. 
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